Non-Gaussianity and Excursion Set Theory: Halo Bias 
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We study the impact of primordial non-Gaussianity generated during inflation on the bias of halos 
using excursion set theory. We recapture the familiar result that the bias scales as k~ 2 on large 
scales for local type non-Gaussianity but explicitly identify the approximations that go into this 
conclusion and the corrections to it. We solve the more complicated problem of non-spherical halos, 
for which the collapse threshold is scale dependent. 



I. INTRODUCTION 

The standard cosmological model that fits a wide va- 
riety of observations is based on inflation, in particular 
on the production of a nearly scale-invariant spectrum 
of adiabatic perturbations at very early times. Inflation 
also predicts, and this too is verified by the data, that 
the perturbations should be nearly Gaussian. While in- 
flation is successful in explaining the current suite of ob- 
servations, it has not been successfully connected to the 
rest of physics. Equivalently, the mechanism that drove 
inflation has not been identified. For these purposes, 
the small deviations from scale-invariance or Gaussianity 
may prove crucial The simplest single field slow- roll 
models of inflation predict negligible non-Gaussianity, so 
a detection has the potential to rule out a large class 
of models. Multiple-field models, on the other hand, of- 
ten predict levels of non-Gaussianity within the range of 
upcoming surveys. 

The cosmic microwave background is the most obvi- 
ous place to search for non-Gaussianity, as the perturba- 
tions are observed when they are still small, and there- 
fore relatively unprocessed. Large scale structure, on 
the other hand, is comprised of highly evolved pertur- 
bations that are nonlinear and hence different Fourier 
modes have mixed with one another. Even if the primor- 
dial perturbations were perfectly Gaussian, the observed 
large scale structure would be highly non-Gaussian. The 
trick to using large scale structure is to identify a feature 
in the spectrum that can be caused by primordial non- 
Gaussianity but not by standard gravitational instabil- 
ity. Recently, such a feature has been identified Q in the 
form of scale-dependent bias0 While it has been known 
for some time that primordial non-Gaussianity affects the 
abundances of rare objects , the signal is expected to 



1 Generically effects that scale as k~ 2 , mimicking the scale de- 
pendent bias on large scales resulting from primordial non- 
Gaussianity, can arise from rclativistic effects in general relativity 
with Gaussian density fluctuations. However, the effects due to 
primordial non-Gaussianity are much stronger than those arising 
dynamically from relativistic effects (see e.g. Udl)- 



be diminished due to the gravitational evolution generi- 
cally evolving the distribution away from Gaussian. 

The initial argument for scale-dependent bias was 
based on both simulations and the statistics of hi gh p eak 
regions pi. [Iol - [l6j . Subsequently, several groups [17l - fl9| 
have applied the peak-background split in the context of 
excursion set theory (20j . Here we apply a recent gen- 
eralization of the excursion set approach [2l|, [22| to the 
problem of the clustering of halos in order to extract the 
scale-dependent bias term and to understand under what 
conditions it holds. We extract the large-scale k~ 2 be- 
havior of the bias for the case of spherical collapse and 
also generalize to ellipsoidal collapse. Here too the bias 
contains a scale-dependent k~ 2 piece, but the coefficient 
differs from that obtained assuming spherical collapse. 

This paper is organized as follows. In SpH we briefly 
review the path integral approach to the excursion set be- 
fore we derive the conditional and unconditional crossing 
rates in ijllll and mV\ respectively, and finally, the halo 
bias in $V] In WII we extract the bias parameters. We 
conclude in i jVIIl In Appendix [XJ we compare our results 
to a case where the probability density can be evaluated 
exactly - where the barrier depends only linearly on time. 



II. PATH INTEGRAL APPROACH TO THE 
EXCURSION SET 

In this section we briefly describe the path integral ap- 
proach to the excursion set and establish our conventions. 
The excursion set theory was developed in a seminal pa- 
per by Bond, Cole, Efstathiou and Kaiser [2(|. For a 
nice, accessible review, see Zentner (23|. The path in- 
tegral approach to the excursion set was developed in a 
series of papers by Maggiore and Riotto [2l|, [22|, [24| and 
extended by De Simone, Maggiore and Riotto to include 
moving barriers and conditional probabilities [111 [2(| . 

As is usual in excursion set theory, we consider the 
density fluctuation with respect to the average density p, 



<J(x) 



- P 



(1) 



smoothed over a region of radius R: 

a H (x) = J dVw fi (x-x')i(x') = y AL^ fl (fe)«5 ke *-. 

(2) 

In this expression, Wr(x — x') is a window function with 
characteristic radius R, and Wn(k) is its Fourier trans- 
form. Translation invariance implies we can choose x, 
and thus for convenience we take x = 0, and suppress 
the argument on the smoothed density field from now 
on, 5ft = Sr(x = 0) . Following the notation of [2l|, we 
denote by Sr = a 1 (R) the variance of the smoothed field 
Sr, 

r°° p ~ 
S R = J d\nk^\W R (k)\ 2 P(k) (3) 

where P{k) is the matter power spectrum. As shown by 
Bond et. al [2(|, if the density fluctuations, Sr, are purely 
Gaussian, and the window function is a top-hat in mo- 
mentum space, Wi?(/c) = 6(fc_R — k) where kR = 1/R, 
then the evolution of Sr. when considered a function of 
'pseudo-time', Sr, is Markovian. That is, each time step 
is uncorrelated with the previous one, and the future evo- 
lution of Sr is independent of its history. In this case one 
can show that the probability density that the field Sr 
takes a value Sr> when smoothed on a scale R' satisfies 
the Fokkcr-Planck equation, with the variance, Sr, play- 
ing the role of pseudo-time (see [23[ for a review). In 
order that one does not suffer from the so-called 'cloud- 
in-cloud' problem, an absorbing barrier condition is im- 
posed as a boundary condition on the probability density 
at the critical threshold for collapse. This means that the 
probability distribution at time Sr includes only trajec- 
tories that never crossed the critical threshold at earliei0 
times Sr> < Sr. 

Obtaining the probability distribution in the Gaussian- 
Markovian limit simply amounts to solving the Fokker- 
Planck equation with various boundary conditions. How- 
ever, in the general non-Markovian limit one must solve 
a complicated i nteg ro-diffcrcntial equation for the prob- 
ability density [2l| . Moreover, even in the Gaussian- 
Markovian limit, it is difficult to obtain solutions for the 
probability distribution in the case where the absorbing 
barrier conditions are more complicated than linear in 
the pseudo-time variable. 

Until recently it was analytically intractable to move 
far beyond the constant collapse threshold and sharp k- 
space filtering assumptions of the original formulation 
of the excursion set [2fJ. Using the original techniques, 
more complicated collapse thresholds have been consid- 
ered in order to describe collapse scenarios beyond spher- 
ical, and the formation of ionized regions during reioniza- 
tion [27], [Hj]. However, in order to capture the effects of 



2 Recall that the variance monotonically increases with decreasing 
R, so "early" times correspond to large scales. 



non-Gaussian fluctuations and more physically realistic 
filter functions, new techniques for solving for the prob- 
ability distribution when the variable Sr does not evolve 
in a Markovian fashion with Sr were needed. In recent 
years there have been several breakthroughs along this 
avenue [HI, M, MM, ■ 

In this work we will follow the path integral approach 
of Maggiore and Riotto [111, [H| ■ Their formalism is gen- 
eral enough to account for non-Gaussian correlations be- 
tween the fluctuations [22^ . for critical thesholds that 
evolve with the smoothing scale and for window func- 
tions WR(k) other than k-space top hats. Rather than 
attempting to solve a differential equation for the proba- 
bility distribution, in the approach of Maggiore and Ri- 
otto one directly constructs this quantity by summing 
over paths that never exceeded the threshold, that is, by 
performing a path integral. 

Following Maggiore and Riotto, we consider the 
smoothed density field, S(S) (= Ss R ), as a stochastic vari- 
able with zero mean (S(S)) = 0. This variable evolves 
stochastically with time S, and we refer to single re- 
alization of a S(S) as a trajectory. We then consider 
an ensemble of trajectories of this stochastic variable all 
starting from the same initial point S(So = 0) = 0, and 
we follow them for some time S. We discretize the time 
interval [0, S] into steps, AS = e, so that Sk = ke, where 
k g {0,1,.., n}, and S n = S. A discretized trajectory 
then, is defined as a set of values {Si,S2, ...,<5n} so that 
S(Si) = 5 t . 

In this path integral approach, the fundamental object 
is the probability density in this space of trajectories, 
which is given by 

W(S ; <5i, S n ; S n ) = (S D (5(S 1 ) -&)... 5 D (5(S n ) - 5 n )) 

(4) 

where S D is the Dirac delta, which we explicitly denote to 
avoid confusion with the density perturbation. Note that 
W is defined for all possible values of S including those 
above the threshold. The bclow-threshold requirement 
is enforced at the next stage when computing Tl(S n , S n ), 
the probability of arriving at point 5 n at time S n through 
trajectories that have never exceeded some threshold, 
B(Si) = B l 

/B]_ pB n -i 
d8 x ... / dS n _ 1 W(S ;5i,...,5 n ;S n ). 
-00 J —00 

(5) 

The upper limits on each S integral allow for the possibil- 
ity of scale-dependent barriers. Although the threshold 
value of S is constant for spherical collapse (equal to 1.686 
today), in non-spherical collapse and other applications, 
the barrier is generally scale-dependent. 

The fundamental quantities, W and n, are the ingre- 
dients needed to connect to observations. The fraction 
of the universe contained in collapsed objects on scales 
larger than S n is found from Eq. by integrating over 



2 



all below-threshold values of 8 n and taking the comple- 
ment: 



/coll(5n) — 1 



dS n U(S ; 5 



(G) 



Differentiating this gives the formation rate of objects 
collapsing on the scale S n : 



dfcoll(Sn) 

dS n ■ 



(7) 



This crossing rate is the unconditional rate that allows for 
trajectories passing through any intermediate values of S. 
It is also useful to compute a conditional rate that fixes 
one of the intermediate values, as F(S n \5 m , S m ), where 
the second arguments dictate the time and value of the 
intermediate fixed point. The conditional crossing rate 
is found from the analogous equation to Eq. ([6]), where 
the unconditional probability that a trajectory reaches 
a point S n is replaced by the corresponding conditional 
probability that a trajectory reached a point S n at time 
S n having had a value S m at intermediate time S m [26], 

n(<5 n , s n \s m , S m ) = 

I-L dSl ■ ■ ■ dSm ■ ■ ■ J-^dfin-iWiSo; <Si, S n ;S n ) 



J-^dSi . . . J ™ l d8 m - 1 W{5o; Si, S m ; S m ) 



(8) 



The dS m in this expression indicates that we do not in- 
tegrate over the variable S m . 

The ratio of the conditional and unconditional rates 
quantifies the fact that halos (on scale S n ) are more likely 
to form in overdense regions (on scale S m < S n ) than in 
an average place in the universe. The halo overdensity 
on large scales S m in initial Lagrangian space is (33| 



1 + 5; 



halo _ J~(S n \S m , S m ) 



(9) 



In this work, we are interested in the halo bias b(k), 
which relates the halo overdensity to the matter overden- 
sity, 



S ha[ °(k) = b(k)S(k). 



(10) 



Extracting b(k) from the smoothed quantities used in 
the excursion set formalism will require a little work, but 
the basic idea is that quantities smoothed on scales R 
carry information about wavenumbcrs k of order R . So 
we will expand the right-hand side of Eq. (|9|) in S m and 
identify the coefficient of the first order term. The haloes 
are collapsed on scales S n , but we are interested in their 
clustering on the much larger scales associated with S m . 
This coefficient can then be massaged to extract b(k). 

As is, the expression for the probability density, 
Eq. (01, is not very useful. One of the insights of Mag- 
giore and Riotto was to manipulate the terms so that 



the right side of that equation contains sums of p-point 
functions. The trick is to use the integral representation 
of the Dirac delta function, 



S D {x) 

so that 

W(S ;Si,...,S n ;S n ) = 
where we have defined 

poo 

v\ = 



d\ 



-iXi 



(11) 



T>\ e l £"=i XiSi (e' 1 ££=i ^(S*)), 

(12) 



d\i d\ n 



(13) 



, (27T) (27T)- 

The expectation value, (e _l ^ =1 A i 5 ( s i)) can be rewritten 
as 

( e -*£?=iA*«(S«)} ( 14 ) 



exp 



E 

p=2 



(-i)P 



J2 ^-^MSji)-ASi r ))c 



where (S(Sji)..S(Sj )) c is the connected p-point function. 
In the Gaussian case, the p-pt functions with p > 2 van- 
ish, and the probability density reduces to the limit: 



W, 



On 



VX e'S-i^^^^i XiXASiSi)^ ( 15 ) 



where the subscript denotes the initial and final times 
of the trajectories; the superscript identifies this as the 
Gaussian limit; and we suppress the dependence on the 
Si, ... , S n -i. We will also have occasion to invoke trajec- 
tories with starting time S m instead of S^; these will be 
drawn from the distribution W / T f i „ and will be integrated 
to obtain the probability 



nf 



n s (5m) S m ; S n ; Sn). 



(16) 



Note that the Gaussian expression in Eq. (fT5|) satisfies 
the identity 

/oo 
VX X kl ... X kn (17) 
-oo 

where di = d/d8i. Using this identity in Eqs. (TT2"j) and 
(jT3J) leads to a general expression for the probability den- 
sity in terms of the p-point functions and derivatives of 
the Gaussian density. 



W 0n = W(S ;Si,...,S n ;S n ) 
(-!)*> 



exp 



E 

p=3 



3U — Jp = 1 



(18) 



3 



The effects of non-Gaussianity primarily arise from 
the non-zero three-point function, so we will drop all 
higher order p-point functions (and all terms quadratic 
and higher in the 3-point function). In this case, the 
probability density reduces to 



TT'c 



Or? 



1 ™ 

1 - - ^2 ($i8jS k } c didjd k 



w§ n . 



(19) 



We are never interested in the raw probability density 
described by Eq. ([TO]) . Rather, we are interested in the 
probability of arriving at a point S n at time S n or in the 
collapse fraction. Thus, rather than evaluate Eq. (fT9|) 
directly, we will integrate over intermediary values of 8i 
to form II and the (observable) functions that can be 
constructed from it. 

One well-known limit is the case when the barrier is 
constant (all i?'s in Eq. (|5|) are the same) and the win- 
dow function in Eq. is a top-hat in fc-space. The top 
hat filter in fc-space means that moving to larger S cor- 
responds to adding more Fourier modes, and since each 
Fourier mode is independent, this filter leads to a Marko- 
vian trajectory. Throughout we stick to the top-hat in 
fc-space filter. In this limit of Gaussian perturbations (g) 
and constant barrier (cb), the probability of reaching 6 n 
at S n starting from S = at S = is 



n s ' ch (s n ,s n ) = , 1 f e- b -i 



-<5^/2S„ _ e -(2B-8 n y /2S r , 



(20) 

The observables then are easily computed: the collapse 
fraction via Eq. (j6|), 



/coll (Sn) 



erfc 



B 



and then the formation rate from Eq. ([7]), 

B 



(21) 



(22) 



For the remainder of this work, we will stick with the top 
hat in fc-space filter, but will work to extend the collapse 
fraction and formation rate to allow for a moving barrier 
and a non-zero 3-point function. In our notation, this 
means we will work to generalize and drop the s and cb 
superscripts!! 



III. UNCONDITIONAL CROSSING RATE 



Starting from Eq. (|19[) . we compute the probability 
of reaching S„ at time S n allowing for any possible set of 



3 Note that the use of a top hat filter in k— space means that 
in the Gaussian limit the trajectories become Markovian, and 
thus the limit we refer to as the Gaussian limit is technically the 
Gaussian-Markovian limit. The effects of smoothing with a more 
realistic filter are relatively straight forward to include |2 ll |2{J , 
see |2S| . |31| . I32I for a study of their effects on the halo bias. 



intermediate <5's. To carry out the calculation beyond the 
Gaussian, constant-barrier limit of Eq. (|20|) , we first need 
to specify how the barrier depends on scale. Following 
[26j], we can specify any barrier by its Taylor expansion 
about the endpoint 



B, — B„ 



£ 

P =i 



B 



(p) 



-(Si S n ) p . 



(23) 



A simple way to account for the moving barrier when 
computing II in Eq. ([5]) is to shift integration variables 



5' 



00 ( P ) 



£ 

P =i 



(S> - s n )p. 



(24) 



Then all the upper limits in the integrals leading to II 
become identical, equal to B n . The integrand (Eq. ([19)) ) 
is an operator acting on Wy n , which now becomes: 



/OO 
V\ cxp z^A 
-°° I 1=1 

1 n 1 

-5 £ XiXj(SiSj) \. 

) 



00 (p) 



6 t + J2^r^-s n r 



p=l 



(25) 



Integrating to form II leads to a Gaussian and non- 
Gaussian term: 

1 - 

m(8 n ,S n )=m(S n ,S n )-- ( S i § j § k)c 



G 

n/ M 

.1=1 j -° a 



i,j,k=l 



didjdkW^. (26) 



A. Gaussian limit 

Consider the Gaussian part of this, the first term on 
the right in Eq. (|2"6"|) . In the limit of a constant barrier, 
this is straightforward to obtain and one finds the result 
Eq. (|20|) . However, for a moving barrier, things are com- 
plicated by the presence of the additional factor in the 
exponent. To evaluate the probability distribution in this 
case, we expand the exponential 



{n 00 (p) 

i £ A *£n-^- 5 «) p 

n 00 R (p) 

= i+££^-w 

i=i P =i 1 



(27) 



n 00 (p) ( q ) 

1 v - ^ -Dn -°n 



2 p \ q \ 



(Si - S n ) P (S 3 - SnYdidj + . . . , 



where the derivatives arc understood to be acting on the 
constant barrier limit of W. In the Shcth-Tormen ap- 
proximation (27l. |34|. one approximates (S n — Si)?^ 1 ~ 



4 



S^" 1 while simultaneously truncating the sum at p = 5 
terms, as discussed in [261 ] . One can then evaluate the 
unconditional probability in the presence of a moving 
barrier in the Gaussian limit. This is given by (25l. |2o| 



IF(6 n ,S n )~IP' cb (6 n ,S n ) 



+ 2 



. Br, 



V^s 3 



5 n ? ^21 
e 



'V 0n 

; — v 2 

1 On 



where 



E 

P =i 



Pi 



(28) 



(29) 



and the . . . refer to terms higher order in an expansion 
in Von, and we are assuming that the barrier is varying 
only slowly in time, S„. It is straightforward to integrate 
this over 8 n to obtain the collapse fraction, 



/ c g oll (S„) =erfc 



B. 



2S T: 



On 



V 2 

1 On / H 2 

Or, 



V 2 

1 On 



-Bl/2S„ 



erfc 



V 25„ 
B n 



'-cvfc 



B„ 



S 2 



S=B„ 



(30) 



To obtain the formation rate, or the mass function, there 
are two ways of proceeding. One could directly differ- 
entiate the expression in Eq. (|30[) with respect to S n 
to obtain the formation rate. Alternatively, since the 
Gaussian-Markovian probability distribution satisfies the 
Fokker-Planck equation, the formation rate can be found 
directly by differentiating the probability distribution, 



(31) 



Of course, these two prescriptions should lead to the same 
result. However, due to the rather crude approximation 
that leads to Eq. (|28|). this expression no longer satisfies 
the Fokker-Planck equation. This is an obvious draw- 
back of the crude approximation we used in order to 
evaluate the probability distribution. In Appendix we 
demonstrate that enforcing the Fokker-Planck equation 
when calculating the formation rates (i.e. using Eq. (|31|) ) , 
rather than directly differentiating the collapse fraction, 
leads to results that are consistent with those obtained 
by an exact treatment of the probability distribution in 
the case of a linear barrier where no such approximation 
needs to be made. Thus, for the rest of this work, where 
possible we apply the Fokker-Planck equation. 



Using Eq. ([3T1) and Eq. (|28|) we obtain the uncondi- 
tional formation rate for a general barrier 



-F g (S„) 



B n + Vor, 



(32) 



Eq. (|32[) is a compact expression that agrees with that 
obtained in previous work, but as described above, it 
hides a subtlety in the Sheth-Tormen approximation. 
Had we simply differentiated Eq. (f30|). we would have 
arrived at a considerably more complicated looking re- 
sult, however, numerically the difference is small. 



B. Non-Gaussian contribution 

To compute the non-Gaussian piece in Eq. (|26p . we 
need to carry out the sums (integrals in the limit that 
e — > 0). Maggiore and Riotto showed that the simple ap- 
proximation of evaluating the 3-point function at the end 
point i = j = k = n leads to the first term in an expan- 
sion in the parameter S n / B 2 . A useful way of evaluating 
the remaining integrals is to use the identity [22|, [2j| 



d 3 " 
f cotton) =~ 

jl J2 J3 = l 



dBl 



dSi 



i=i 



(33) 



This is almost identical to the sum and integrals in the 
last term in Eq. (|26[) . with the exception that the right 
hand side in Eq. (|33[) also includes an integral over S n . 
To use the identity then, we can integrate Eq. (|26p over 
6 n (which we need to do anyway in order to compute the 
collapse fraction). This leaves 



ng 1 3 ® 3 

/coll = g^™' Qg3 /colU^n) ' 



(34) 



To obtain the crossing rate, we differentiate this with 
respect to S n . The derivative acting on the three point 
function is straightforward. To evaluate the derivative 
acting on /^ oll , we pull the derivative all the way through: 
/ dSdU(S; S n )/dS n ; use the Fokker-Planck equation and 
then integrate by parts so that the non-Gaussian part of 
the collapse fraction is 



^(Sn 



We have defined 



1 - tS% 



d 3 \ dm 
dBj) ~W 

if' 



5=B„ 



(2S n S 3 + S 2 S' 3 )—fUS n ). (35) 



S, 



•5 n 



(36) 



and denoted derivatives with respect to S n by primes, '. 
Using Eq. $5U§ and Eq. (J2SJ we find that the full uncon- 
ditional crossing rate to leading order in the three point 



5 



function is given by 



, 5 3 e ^ 



(37) 



B n - &PonSn + 75 r 2 j 



'PonBn (-B„ + S n ) — lOB^Sn 



3y/2nS„ 



[-^n — B n Von — S„ + 2V Qn \ 




B n -S m 



IV. CONDITIONAL CROSSING RATE 



SL -SL 



In order to evaluate the halo ovcrdensity, Eq. ©, we 
also need the conditional crossing rate. The major dif- 
ference between evaluating the probability distribution 
in this case (Eq. (|5J) and the previous case (Eq. (j^}) is 
that the intermediate point on the trajectory denoted by 
S m is not integrated over. Using the relation between W 
and W g in Eq. ([TO)), the starting point for the conditional 
crossing rate calculation is then 



I Sm : S m ) 



i,j,k=l 
Bi 

d5i . 



d5\ . . . dS m . . . dS n -i 

i J - 

1 \ 

1 - r E (SiWkUdidjdk w§ n 



i,j',fe=l 



Gaussian contribution 



FIG. 1. Top panel shows a trajectory from large scale S m 
to small scale 5*„. In the Gaussian, Markovian limit, the 
probability for this trajectory does not depend on prior steps, 
so is equivalent to the trajectory depicted in the bottom panel 
with the shifts enumerated in Eq. (140 \ . 



evaluated in Eq. (|28|) ) as long as we make the substitu- 
tions 



(40) 



For example, the first term in Eq. (f28|) - the constant 
barrier piece - becomes 



S n - 


■> 5 n 


- s m 


B n - 


■> B n 


- Sm 


S n ~ 


^s n 


~ S m 



n Si {S n i S n | S m , S m ) 



x ( e -[<5„-5 m ] 2 /2[S„-S m ] 



1 



^/27r[S' n — <S m ] 

-(2B~S n -5 m ) 2 /2[S n -S m h 



(41) 



We first consider the Gaussian piece of Eq. ([55} by 
dropping all terms with 3-point functions in them. In 
this limit, there is a simplification due to the Markovian 
nature of the trajectories. Namely, the probabilities do 
not depend on the paths chosen but simply on the end 
points. So Wq u in the numerator can be replaced with 
W ( f m W / ^„. Then all the integrals over S 1 ,S 2 ,---, 5 m -i 
are identical in the numerator and denominator, leading 
to 

/B m +l fSn-i 
dSm+1 ■■■ dS n -lW mn . 
-oo J — oo 

(39) 

As depicted in Fig. [TJ this expression is identical to the 
unconditional Gaussian Markovian probability (which we 



B. Non-Gaussian contribution 

To evaluate the non-Gaussian corrections to the condi- 
tional collapse rate, we again rely on the approximation 
of considering the three-point functions only at their end 
points. In the case of unconditional collapse, this trans- 
lated to setting (SiSjSk) — > (S n S n S n ), thereby simplifying 
the sum over i, j, k. Here, too, we set all three-point func- 
tions to their values at the end point but when i < m, 
the end point is time S m , not S n - For example, 

n m— 1 n m— 1 

E EvW*>c«<<wOc E E- ( 42 ) 

i,j= m+1 k—1 i.k— m+1 k—1 
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and working to linear order in the long wavelength mode, Rotating the integrals in the analogous manner as Eqs. 

(|23)) -([24 |) and making use of identities such as Eq. (|33|) 
as well as 



5 m , we can write (see also |35|) 

n 



m— 1 n 



fc=l *,j=m+l 
n 

n 



m-1 ,,5 



dB 



N 



n§ m = E 



(44) 



(43) 

we find that the non-Gaussian part of the conditional 
collapse fraction can be written 



/coll {£>n\fimi S m ) — 



\°n)c dB a 



3(^(5,, 



d_ 

dBldL, 



( 9 







\ dB n dS r 



Inn? 



'dBl 



/coll ( &n 1 &m j <S m ) , 

(45) 



where / s j^SVil^m, S'm) is the conditional collapse fraction in the Gaussian limit, which follows from the results of 

Taking the derivatives and working in the limit S m — > with B n S m , we find that the non-Gaussian contribution 
to the conditional collapse fraction is 



/coil (S n \5 m , S„ 



dm 



1 

2 exp 



g 2(S„-S m ) 



(-B„— 
2(S„-S m ) 



^2^(5„ - S m ) 3 

\P'mn ( )) 



(B n -5 m ) 2 { (B n -6 m ) 
(S n S m ) {S n S m ) 



rfc 



B„ — S„ 



\/2(S'n — S m ) 



V 2 



(S n S m ) 



The key result here is the appearance of the combination 

d 



(S m S 2 n ) c (J- 



d6 r , 



111 ii 0.m 



(46) 



(47) 



which multiplies the last line in Eq. (|4"6"|) . The appearance of the terms proportional to ((5 TO <5 2 ) C in Eq. (|4"6"]l encode the 
fact that the conditional probability of collapse now depends on the details of the correlation between long and short 
wavelength fluctuations. We will see that these correlations give rise to a scale dependent bias. In order to extract 
the bias we will eventually correlate the expression in Eq. (|46[) (to linear order in 5 m ) with the underlying density 
fluctuations smoothed on the scale R m , before taking the limit that S m — > 0. While the quantity in Eq. (|47|) remains 
finite after this procedure, the terms in the first line of Eq. (|46|) proportional to (<5 m <5 2 ) c vanish, and we thus drop 
them in what follows. 

We can then evaluate the conditional crossing rate, which is the conditional rate at which objects form at time 
S n , given that they had an overdensity S m at time S m . As we have stressed above, we do not proceed by simply 
differentiating Eq. (|4"6")l with respect to S n . Rather, we write the non-Gaussian contribution to the conditional crossing 
rate as 



mb.NG 



(S n \5 m> S m ) 



^/coll (Sn\8m, Sm ) 



dS n 



d 3 

dBl 



3 y S m S n S n 2 r 



8 d \ d 2 

dB n + dd m r ll °> m dBi 



{2S n Sz + S n S' 3 ) 



f . d 3 



dBl 



-Z\J ' S m (S n 2 m + S n S' n 2 r 



1 dU s 

2 dS n 

d 

dB n dS m 







inn 



0,m 



dBl 



f! n(Sn\S m ,S m ), (48) 
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where we have made use of the Fokker-Planck equation in the second line and primes denote derivatives with respect 
to S n . We have dropped the terms which vanish in the S m — > limit as described above and defined 



(49) 



We then evaluate to find the conditional crossing rate 



J~ m ' (S n \&rm S m ) — — —S' n ; 



g 2(S„-S TO ) 

y/2n(S n — S m ) 



S„3 e 2(S„-S,„) 



[(S n - S m ) - (B„ - 5 m f + (B„ - 5 m )V nm - 2V 2 m ] 



r Pmn{B n — S m ) ((B n — 5 m ) 2 + (S„ — S m )) — 10(-B„ — S n ) 2 (S„ — S m ) 



6 ^{S n - S m ) 5 

+ (B n — 8 rn ) A — 8V mn (S n — S m ) + 7(S n — S m ) 2 



SmS^: 



7m<-> n 2 m „ 
On 



V S m S n 2 7n 



g 2(S„-S m ) 

^/2-rr(S n — S m ) 

(B„-6 m ) 2 
g 2(S„-S m ) 



s m yj2it{s n - s m y 



rfc 



ij^mn (B n S m )) 



B n -8„ 



y/2(S n -S m ) 



(S n S m ) 



{B n — s m y 

Sn Sm 



(B n — S m + ^Pmn) Q'Pmn (Bn Sm) 



2ir {B n -s m ) 2 



S n S, 



Bn &rn 
\/2(S n S m ) 



(50) 



V. HALO BIAS 



We now have all the ingredients with which to extract the halo bias. Taking the ratio of Eq. (|50|) and Eq. (|37|) . and 
expanding to linear order in S m , we find that the relationship between the halo overdensity and the matter overdensity 
is 



S n B n + Vqt 

1 S„3 

6 (B n + Von) 2 

SnS' n3 I 
i(B n + Von) 2 I 



20V Qn B n + 10B 2 n - 9P 2 „ - 3^- (B n + V 0n ) 2 + 7S n \ S„ 



(B n - Von)(B n + 3V 0n ) + S n S r , 



. i B n 



S n \/ S m S' n 2 m 1 



B n 1 A 3Von 

S n B n + Von J B n + Po 

2V 0n 



2 * ^ erfc 



B n + Von V S n B n + Vq 



S n B n + Vo r, 

B„ 



B„ 



'2S n 



2S n 



Srr 
Sn 



2Sn 



(51) 



We can immediately read off from the first line here the 
well-known (27j Gaussian bias in initial Lagrangian space: 



S n B n + Von 

The coefficients of S m in the second and third lines do not 
depend on the large scale m and so represent the scale in- 
dependent contribution to the bias due to non-Gaussian 



correlations between the initial fluctuations. The result 
presented here is slightly different from that reported by 
[2a] who do not employ the Fokker-Planck equation. The 
last two lines contain terms wherein the coefficients of 5 m 
depend on the scale m, so the bias from these terms will 
be scale-dependent. Note that they are proportional to 
the three-point function so vanish in the Gaussian limit. 
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Also, notice that the effects of the moving barrier are 
encoded in the Toon's here. These terms then give the 
full contribution to the scale-dependent bias from non- 
Gaussianity accounting for ellipsoidal collapse. Care is 
required to extract the exact form of the scale depen- 
dent bias, because the bias is typically defined in Fourier 
space, while we have worked with real space fluctuations 
smoothed over various scales. 



VI. THE SCALE DEPENDENT BIAS FROM 
EXCURSION SETS 

We can now extract the scale dependent bias from our 
excursion set result. For simplicity, we begin with the 
case of spherical collapse where the barrier is constant, 
B(S n ) = 5 C . Focusing on the scale dependent parts of Eq. 
(|51[) . and assuming a spherical barrier so that Va n = 0, 
we have 



rhalo /q q <v u i 



The problem is to transform this relation, between the 
real space halo overdensity smoothed on a large scale S m 
with the real space matter overdensity smoothed on the 
same scale, into a Fourier space relationship: 




S h&lo (k) = b(k)5{k). 



(54) 



In terms of the Fourier space fields, the relationship 
between the halo field and the matter field, as derived in 
the excursion set, Eq. (fS"3"f becomes 



d 3 k 

(27T) 5 



S(k)b SD ' sph (k)W(k,R ri 



Srn S n 2 



5 1 



1 



d 3 k' 5(k')W(k,R m ) 

Wf s m 



(55) 



where the SD ' s P h superscript reminds us that we are deal- 
ing with the scale-dependent part of the bias only here 
and, for now, focusing on the spherical collapse limit. 
Now correlate both sides of Eq. (|5"5"|) with 5 m . The left 
hand side becomes 

km 

" P s (k)b SD ^ h (k) -- 



(27T) 



d\nk b SD ' sph (k)A 2 s (k), 
(56) 

where A 2 (fc) = fc 3 P(fc)/27r 2 , and we have defined the 
power spectrum in the usual way, 



(S(k)S(k')) =P 5 (k)(2iry5 d {k-k' 



(57) 



The right-hand side becomes the set of coefficients in 
parenthesis in Eq. (|55|) multiplying 



1 



d 3 k 

(2tt) 3 



Ps{k) 



1 



(58) 



as the k integral here is equal to S m for our chosen filter. 
So the remaining equality is 



dlnkb SB ' sph (k)A 2 s (k) = 




To proceed we need information about the shape of 
the three-point function that enters into Eq. ([59)) . The 
primordial 3-point functions are most easily written in 
terms of the gravitational potential, so first recall that 



where 



5(k,z) = M(k,z)^{k) 
2 k 2 T(k)g(z) 



M(k,z 



3n m H 2 (l + z) 



(60) 



(61) 



and here T{k) is the matter transfer function normalized 
at unity as k — ► 0, g(z) is the linear growth of the grav- 
itational potential during the matter dominated epoch, 
and $(k) denotes the primordial value of the gravita- 
tional potential (hence no z dependence) . A generic 3-pt 
function of $ is characterized by the bispectrum, 

($(ki)$(k 2 )$(k 3 )) = B 9 (ki, fe, fc 3 )(2^) 3 5 3 (ki + k 2 + k 3 ) 

(62) 

We can then evaluate the 3-pt function of density fluctu- 
ations, 



1 dki f kn dk 2 



dcos 6 



(63) 



x ^^M(k 1 )M(k 2 )M{k 3 )B 9 (k 1 ,k 2 ,k 3 ), 

where 9 is the angle between k\ and k 2 and we have as- 
sumed that \ks\ = \ki + k 2 \ < k n . Following the notation 
of [H, dl , w e introduce the function (not to be confused 
with the crossing rate!) 



1 



S n A%(k) 



k 2 dk 2 
2tt 2 



M(k 2 )M(k 2 -k\) 



x B 9 (ki,k 2 ,\k 2 -k\), (64) 



where k = kz, and z, is an arbitrary unit vector. In the 
limit k m <C k n , Eq. ([59]) then reduces to 



A 2 (fc) 



d\nkA 2 s (k)b SB ' sph (k) = dink 



M{k) 

(65) 



At this stage, we have an equality between two sums 
of things. To see that we can equate each term in the 



9 



sum, imagine that we measure these sums for two slightly 
different scales, k m and k m + Sh, then we can extract 
information about the contribution to the sum from scale 
k m by differencing these quantities. This then implies, as 
well as being an equality between sums, equality between 
each term in the sum. Then the scale dependent bias can 
be read off: 



& SD > s P h (fc) =M-\k) 



5 1 



1 Jt{k) 



d In S„ 



(66) 



which agrees with the result derived by [18|. Note that 
we naturally obtain the important second term, which 
can be thought of arising due to the fact that a scale 
dependent rescaling of the variance also changes the sig- 
nificance interval that corresponds to a fixed mass. As we 
will see, this extra term vanishes for the local ansatz, but 
is non-zero and generally important for other bispectrum 
shapes [HI, HH and in particular largely ameliorates the 
discrepancies noted by [37| . 
For the local ansatz, 



so that 



K « 4/, 



XL 



(68) 



a constant. Therefore only the first term in Eq. ([66 
contributes leaving 



^SD, sph, local 



2./': 



NL 



M(k) 



5 1 

Sri 



(69) 



for the scale dependent bias in the spherical collapse case, 
which agrees exactly with the result of 0]. 

Using the same techniques, we can extract the scale 
dependent bias from Eq. (|BT|) for the full ellipsoidal col- 
lapse problem. It too has pieces proportional to both J® 
and dJ-^/dSn- We separate these out as: 



b SD (k) 



1 dJ*(k) 



2M(k) M(k) dlnSr, 



(70) 



k 2 ,k 3 ) = 2/ NL (P(fc 1 )P(fc 2 ) + perm.) (67) where the coefficients are: 



and 



■"a — 



d n = 1 



B n + Vqt, 



ZPor. 



Fhr V 2 J*L 



V 2 B l 

-e^crfc 




(B n + Von) V S n (B n + Von) 



B,, 



B,, 



(71) 



(72) 



These equations are the main results of the paper. 

Ignoring the terms proportional to dJ-f^k) / dS n in 
Eq. ([70)1 (which vanishes for the local ansatz) we see that 
the coefficient of the k~ 2 term, c„, has a piece equal to 
the barrier times the Gaussian bias (B n b s ). This is the 
standard coefficient. However, the full treatment here 
has uncovered other terms, the remaining ones on the 
right in Eq. (|TTj) . In Fig. [5] we plot c n for the spherical 
collapse barrier (i.e. B(S n ) = S c ) and for the ellipsoidal 
collapse barrier of [38^ : 

B(S n ) = A[l + /^Hi (73) 

where v = 5 2 /S n , a w 0.7, a w 0.6 and f3 rj 0.4. For the 
spherical collapse barrier, c n is given by the well known 
expression c„ = S^/S n — 1 (e.g. Eq. (|69|) ). For the ellip- 
soidal collapse model, we must use the more complicated 
expression we have derived in Eq. (|71[) . 

The naive prediction for c„ (which holds in the case 
of the spherical collapse barrier) is that it is related to 
the scale independent bias via c„ = b e B n (S n ). In fact, 
one way to try to extract constraints on /nl is to fit 
the power spectrum with a bias of the form 1 + b s + 



2fisiLb e S c M~ 1 (k). We have shown, however, in Eq. [70] 
that this simple relation does not appear to be correct 
when the barrier evolves with the smoothing scale. In 
Fig. [3] we show the ratio of c„ to the naive prediction 
as a function of S n for the ellipsoidal collapse barrier 
given in Eq. (|75|) . Note that our result depends heavily 
on the approximation of evaluating the 3-point functions 
at their end-point values, as in for example Eq. (|42p . 
This is known to break down as S n becomes large, of 
order B^. Therefore, the upturn in Fig. [3] may be an 
artifact of our approximation and almost certainly does 
not persist as S n increases. Techniques that move beyond 
this approximation [39[ are clearly called for. 



VII. CONCLUSIONS 

Making use of the path integral approach to excur- 
sion set theory, we have derived the bias of collapsed 
objects. We have kept all terms linear in the long wave- 
length density fluctuation and we have allowed for a com- 
pletely general moving barrier. The bias contains both 
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FIG. 2. The coefficient (as defined in Eq. (f70|)) of the scale 
dependent bias, c„, as a function of the smoothing scale S n - 
Solid (black) curve and dashed (blue) curve show the result 
in Eq. [70] for the ellipsoidal collapse model of (i.e. Eq. 
I73|l and the spherical collapse model (i.e. when B n = <5 C and 
Von = 0) respectively. The upper x-axis shows the mass corre- 
sponding to a particular S n assuming a flat, ACDM cosmology 
with n b h 2 = 0.02, n c h 2 = 0.11, h = 0.7 and a 8 = 0.87. 

the scale independent contributions, analogous to those 
first reported by [26j for a general barrier, and contri- 
butions that depend on the long wavelength smoothing 
scale. We have demonstrated how this dependence can 
be interpreted as a scale dependent biasing and shown 
that our result reduces to the famous Dalai et. al. [|[ re- 
sult in the limit of spherical collapse and the local ansatz. 
In the spherical collapse limit for a general bispectrum 
shape, we reproduce the results of Desjacques, Jeong and 
Schmidt [HI, [36| , including the previously overlooked ad- 
ditional term. 

In the case where the barrier is not constant and de- 
pends on the smoothing scale, we find that the coefficient 
of the scale dependent bias is no longer simply related to 
the Gaussian bias parameter, but rather contains addi- 
tional terms (Eq. (fTTj) ) that might affect the extraction 
of /nl from upcoming surveys. While the simple relation 
is recovered on sufficiently large mass scales, on smaller 
mass scales, we find a significant departure from the ex- 
pected result. 

In arriving at this result, we have made several approx- 
imations. Following standard techniques in evaluating 
the effect of a moving barrier (which characterizes ellip- 
soidal collapse), we truncated the probability distribution 
in Eq. (|28[) . This approximation causes II to no longer 
satisfy the Fokker-Planck equation. However, we have 



FIG. 3. The coefficient of the scale dependent bias, c n , from 
Eq. (|70[) for the ellipsoidal collapse model of [3S|], normalized 
to the naive prediction c n = b K B n . The upper x-axis shows 
the mass corresponding to a particular S n assuming the same 
cosmology as in Fig. [21 



shown in Appendix [A] that when used in combination 
with the Fokker-Planck equation to calculate the halo 
bias, this approximation leads to results consistent with 
an exact treatment (at least to cubic order in derivatives) 
when compared to the case of a linear barrier where the 
probability distribution can be computed exactly. For 
this reason, we do not think this approximation leads to 
the upturn in c„ at large S n depicted in Fig. [3J 

A second approximation treating non-Gaussianity in 
the excursion set is likely to have more effect on our fi- 
nal answer. We evaluated the 3-point function at the 
end point of the trajectories, as in Eq. (f4"2j). As noted 
there, one may think of this as keeping the zeroth or- 
der term in a Taylor expansion about the final time S n . 
It is straightforward to calculate the contributions from 
additional terms in this series by making use of the re- 
sults and methods of [2l|, [22|, UM ■ This approximation 
can be shown to correspond to an expansion in Sn/B%. 
Furthermore, using the saddle point techniques of [40( it 
is possible to resum a number of higher order corrections 
to this formula to obtain a more accurate formula for the 
bias. However, to reproduce the existing results in the 
literature, the zeroth order approximation used here ap- 
pears sufficient. We leave the calculation of higher order 
corrections to future work. 

While this paper was in preparation, we became aware 
of the work j39|, which also considers the halo bias in 
the path integral formulation of the excursion set. While 
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the authors of [39| consider only a constant barrier, they 
include next-to-leading order corrections from relaxing 
the approximation at Eq. (J42J . The results of this work 
are consistent with those found in j39| the limit where 
the barrier is constant, and one restricts to the leading 
order result. 



ACKNOWLEDGMENTS 

We thank Neal Dalai, Anson D'Aloisio, Wayne Hu, 
Donghui Jeong, and Toni Riotto for useful discussions 
and the authors of [39[ for sharing an earlier draft of their 
work with us. This work was supported in part by Na- 
tional Science Foundation under Grant AST-090872, the 
Kavli Institute for Cosmological Physics at the Univer- 
sity of Chicago through grants NSF PHY-0 114422 and 
NSF PHY-0551142 and an endowment from the Kavli 
Foundation and its founder Fred Kavli. SD is supported 
by the U.S. Department of Energy, including grant DE- 
FG02-95ER40896. AL was supported in part by the NSF 
through grant AST-1109156. 



Appendix A: The linear barrier 

In this Appendix we compare our results to a case 
where the probability distribution can be evaluated ex- 
actly. Aside from the limit where the absorbing barrier is 
independent of the smoothing scale, the only other known 
exact solution to the Fokker-Planck equation with an ab- 



sorbing barrier condition is in the case where the barrier 
has a linear dependence on the smoothing scale, 



B(S) =B a + B 1 S 



(Al) 



where Bq and B\ are constants. For this barrier, it is an 
elementary exercise to obtain the probability distribution 



U LB (S n ,S n ) 



1 



\/2ttS„ 



exp 







^exp 


2 On _ 



(2B n - 8 n f 
2S n 



2B' n (B n — S n ) 



which one can verify solves the Fokker-Planck equation, 
with the linear boundary condition Eq. (|A1[) . The super- 
script LB here and in the rest of this appendix denotes 
expressions that are only valid for the barrier in Eq. (|A1|) . 

It is straightforward to obtain the Gaussian formation 
rate, 



B n -S n B' n c _^l 
v / 2^53 



(A2) 



which agrees exactly with the result in Eq. ([3"TT) since for 
the barrier in Eq. (|A1[) 

Von = S n B' n . (A3) 

The conditional crossing rate is found as described above 
by simply translating the Gaussian result using the shift 
of variable in Eq. (|40|). We can also compute the halo 
bias analogous to Eq. (fijTj) , 



ehalo.LB 



Bn _ 1 

Sn [Bn S n B n 

1 / 



Bl 



_ 20S n B' n B n + 10B* - %S n B' n f - 3-^ (B n - S n B' n ) 2 + 7S n 
6 [B n - S n B' n y \ s n 



B n — 2S n B' n 



2Sn 



- 8y/toS* /2 B' n 3 e ""'a" " ((B n - 2S n B' n ) (B n - S n B' n ) - S n ) erfc 

S S' ( 

- 3^ -SB' ) 2 I ^ Bn + SnB ' n ^ Bn ~ ^ SnB 'n) + s n- ^S 2 B' 3 {B n - S n B' n ) 
+ 2V2^S 3 J 2 B' 7 fe {B ^: B '" r ((B n - S n B' n ) (B n - 2S n B' n ) - S n ) erfc ( Bn I 



2S n 



c / pi 

^n 2 m f ^n 



B u 



3S n B n 



2 I S n {B n — S n B' n ) (B n — S n B' n 

SnS n 2 m 1 



2V2^Sj 



(B n — S n B' n ) 



-erfc 



B n 2S n B' n 

\/2S n 



(B n ~2S n B / ) 2 



(B n — S n B' n ) 



(B n — S n B' n ) 



■erfc 



B n — 2S n B' n 



<2Sn 



(A4) 
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from which we can obtain the scale dependent part of the bias 



lLB 



(k) 



~2M(k) 
1 



Br, 



S n B' n 



(B„-2S„B^) 2 



dJ*(k) 



SnB' n 



2S n B n 



B„ 



S n B' n 



(B„ -2Sn,B' n ) 2 



M{k) d\nS n 



(B n — S n B' n ) 



B n - S n B' n 



B r , 



-crfc 



S n B' n 



-crfc 



B n 



2S n B' n 



2Sn 



Br, 



2S n B' n 



2S n 



(A5) 



and thus we find exact agreement with Eq. (|51|) and 
Eq. ([70)) up to terms cubic in B' n , the order to which 
we have evaluated the approximation in Eq. (|28|) . We 
thus conclude that the deviation we find at low masses is 



not a consequence of the approximation of the probability 
distribution, but rather represents a break down of the 
approximations at Eq. (|42[) or possibly of the excursion 
set itself. 
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